Astrophysical S-factor of 14 N(p, 7 ) 15 cE 
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We report on a new measurement of the 14 N(p, 7) ls O capture cross section at E p = 140 to 400 keV 
using the 400 kV LUNA accelerator facility at the Laboratori Nazionali del Gran Sasso (LNGS). 
The uncertainties have been reduced with respect to previous measurements and their analysis. We 
have analyzed the data using the R-matrix method and we find that the ground state transition 
accounts for about 15% of the total S-factor. The main contribution to the S-factor is given by the 
transition to the 6.79 MeV state. We find a total S(0) = 1.7 ± 0.2 keV b, in agreement with recent 
extrapolations. The result has important consequences for the solar neutrino spectrum as well as 
for the age of globular clusters. 

PACS numbers: 24.30.-v, 24.50+g, 26.20.+f 



The capture reaction 14 N(p,7) 15 (Q = 7297 keV) is 
the slowest process in the hydrogen burning CNO cycle 
and thus of high astrophysical interest. This reaction 
plays a role of setting the energy production and neu- 
trino spectrum of the sun [2j as well as the age determi- 
nation of globular clusters j3j. Below 2 MeV, several 15 O 
states contribute to the 14 N(p, 7) 15 cross section (fig. 
HJ): a J" = 3/2+ subthreshold state at E R = -507 keV 
(E x = 6.79 MeV), and 3 resonant states: r = 1/2+ 
at Er = 259 keV, and 3/2+ at E R = 987 keV and 
Er = 2187 keV. The reaction was previously studied 
over a wide range of energies, i.e. E cm = E = 240 to 
3300 keV (0 and references therein). The non resonant 
capture to excited states in 15 O led the authors of ref. 
to an extrapolated astrophysical S-factor at zero en- 
ergy of S es (0) = 1.65 keV b. The data for the capture 
into the 15 O ground state were analyzed using the Breit- 
Wigner formalism and indicated an important influence 
of a subthreshold state at Er = —507 keV, leading to 
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S gs {Q) — 1.55 ±0.34 keV b with a deduced gamma width 
of T 7 = 6.3 eV, thus S tot (0) = 3.20 ± 0.54 keV b. 

A reanalysis [j| of the capture data to the ground 
state 4] using an R-matrix approach indicated a neg- 
ligible contribution of the subthreshold state to the to- 
tal S(0)-factor, mainly due to a significantly smaller T 7 
of this state. This indication was supported by a life- 
time measurement of the Er = —507 keV subthresh- 
old state via the Doppler-shift method @ leading to 
T 7 = 0.411q' 3 3 eV, and a measurement via the Coulomb 
excitation method Q resulted in T 7 = 0.95±°;g5 eV. In 
view of the uncertainties in S gs (0) and thus St o t(0), a 
new measurement of the capture process into the 15 O 
ground state was highly desirable, extending possibly the 
low energy limit below that of previous work i.e., below 
E = 240 keV. We report on the results of such measure- 
ments using the 400 kV underground accelerator facility 
LUNA at Gran Sasso, Italy. 

The 400 kV LUNA facility and the setup have been 
described elsewhere 8] . Briefly, the accelerator provided 
a proton current on target of up to 500 /iA. The abso- 
lute energy is known with an accuracy of 0.3 keV and the 
energy spread and the long-term energy stability were ob- 
served to be 100 eV and 5 eV/h, respectively. Near the 
target, the beam passed a liquid-nitrogen cooled shroud 
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FIG. 1: Level structure of 15 0. Above the 14 N(p, 7 ) 15 
threshold only relevant states are shown. 



(to minimize carbon-buildup on target) and an electri- 
cally insulated collimator with a negative voltage of 300 V 
(to suppress the effects of secondary electrons). The 
water-cooled target was oriented with its normal at 55° to 
the beam direction. The target consisted of a TiN layer 
(with a typical thickness of 80 keV) reactively sputtered 
on a 0.2 mm thick Ta backing. The target quality was 
checked frequently at the En = 259 keV resonance: no 
significant deterioration was observed after a bombard- 
ing time of several days. Typically, a new TiN target was 
used after a running time of 1 week. The stoichiometry of 
the TiN layer was verified via Rutherford Backscattering 
Spectrometry using a 2.0 MeV 4 He + beam, resulting in 
Ti/N= 1/(1.08 ±0.05). 

To deduce cross section values in the non-resonant en- 
ergy region, one can use either thin targets (say, a few 
keV) or thick targets (say, a few 10 keV). Due to sput- 
tering effects at low energies thin targets deteriorate sig- 
nificantly in a short time changing correspondingly the 
observed 7-ray yields. In the case of thick targets, the 
sputtering effect on the observed yields is negligible, as 
long as the thickness stays large enough, which was ver- 
ified experimentally (see above). For the chosen thick 
targets, we made an analysis of the line shape of the pri- 
mary 7-rays as described below to deduce a capture cross 
section. 

For the measurement of excitation functions the cap- 
ture 7-rays were observed with one Ge detector (126% 
efficiency) placed at 55° in close geometry to the target 
(d = 1.53 cm was the distance between the target and 
the front face of the detector). In one experiment the dis- 
tance was increased up to 20.5 cm for the determination 
of the detector efficiency and summing effects. In another 
experiment, Ge detectors were placed at 0° (126 %), 90° 
(120 %), and 125° (108 %), relative to the beam axis 
(d = 7 cm) for the measurement of Doppler shifts, exci- 



tation energies, and angular distributions. Sample spec- 
tra are shown in fig. [21 where the influence of different 
background sources can be observed: (a) beam induced 
background, which is mainly due to proton capture on 
n B, 18 O, and 23 Na while at E p = 140 keV (b) the reman- 
ing beam induced background from n B contaminations 
is equal to the cosmic background. 

In a similar way as described in the line shape of 
the 70-capture transition in 14 N(p, 7 ) 15 may be used 
for various purposes. Fig. [3] shows the 70-ray line shape 
observed at Q 1 — 55° and E — 243 keV, i.e. at the low- 
energy tail of the Er = 259 keV resonance; the dispersion 
in the spectrum was 1 keV/channel. The line shape can 
be interpreted as the sum of the resonant and non res- 
onant contributions; when convoluted with the detector 
resolution, the solid curve through the data points is ob- 
tained. In this way, the drop in the 70-ray yield towards 
lower energies reflects directly the drop of the cross sec- 
tion. For the calculated line shape the dependence on 
energy of the 70-efficiency and of the stopping power of 
protons in TiN was included. The high-energy edge 
of the peak contains the information on the incident beam 
energy, hence possible C-build up on target could be cor- 
rected for; this correction was never larger than 2 keV. 
The method to extract the cross section from the 70-ray 
line shape requires further investigation of the following 
points: (i) energy, width, branching ratios, and strength 
of the Er = 259 keV resonance, (ii) energy calibration 
and efficiency of the Ge detector, (iii) summing effects 
in close geometry, (iv) Doppler shift, and (v) effects of 
angular distributions. In the same way the transition to 
the 6.79 MeV state was analysed. Below E = 170 keV 
the information of the data from the secondary transi- 
tion was used in the analysis. For this energy region a 
constant S6.79 was assumed over the target thickness. A 
detailed description of all these points is given elsewhere 
[Til fT3 | . where a complete analysis of the data to all final 
states will be presented. 

The resonance energy could in principle be determined 
from the 70-energy after correction for Doppler shift and 
recoil. One requirement is the precise energy calibra- 
tion of the Ge detector. We have used the information 
obtained from the experiment with the 3 Ge detectors 
(0 7 = 0°, 90°, and 125°). From @ it is known, that 
the E x = 5181 keV state in 15 shows an attenuated 
Doppler shift and the 6172 and 6791 keV states have 
lifetimes resulting in nearly full Doppler shifts. Using 
the known energy of the resonance from the accelerator 
calibration together with Q = 7296.8 ± 0.5 keV [3 and 
calibration points from radioactive sources we have per- 
formed a combined \ 2 fit of the data obtained at the 
three angles. It was necessary to vary also the excitation 
energies of the first three excited states. In addition, we 
varied the attenuation coefficients for the Doppler shifts. 
The Doppler shift data for the E x — 6791 keV state are 
shown in Fig. 0|and all results are summarized in Table 
U We confirm the results of [fj for the first excited state 
in 15 O but could not extract a lifetime for the 6791 keV 
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FIG. 2: 7-spectrum for 14 N(p, 7 ) 15 obtained (a) at E p = 200 keV (Q = 85 C, t = 62 h) and (b) at E p = 140 keV (Q = 222 C, 
t = 396 h). For case (a) beam induced background originates primarily from traces of 11 B, ls O, and 23 Na in the target, while 
for case (b) the beam induced background was of the same order than the remaining cosmic background which is about 7 
counts in the region of the peak for the ground state transition where a total of 50 counts are observed. 



TABLE I: Excitation energies in 15 O, Doppler shift measurements and 7-branching ratios for the 259 keV resonance (E x = 
7556 keV). 
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C N implanted in Ta 
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state due to a nearly full Doppler shift. For all the sub- 
sequent work we calibrated the 7-ray spectra using the 
new excitation energies given in Table [I] 

The detector efficiency was determined using cali- 
brated radioactive sources and the cascade condition for 
the transitions to the first three excited states at the 
Er = 259 keV resonance, as for all of them no other 
decay than that to the ground state was observed. This 
procedure was performed with the Ge detector placed at 
1.53, 5.5, 10.5, and 20.5 cm distances from the target 
in order to determine the summing-in contribution to 



the ground state transition and the summing-out for the 
transitions to the excited states. It turned out, that the 
summing-in yield was about 3.5 times higher than the ac- 
tual ground state intensity at the 1.53 cm position. This 
1.53 cm position was used for the entire cross section mea- 
surements. The efficiency curve was also calculated us- 
ing the GEANT routine [lj] and found to be in excellent 
agreement with observation [Tl| . Branching and strength 
values of the Er — 259 keV resonance were determined 
with low beam current to avoid dead time effects and at 
far distance (20.5 cm) to minimize summing effects from 
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TABLE II: R-matrix parameters for capture transition to the ground state 
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FIG. 3: Typical 70-ray line shape of 14 N(p,7) ls O obtained at 
E p = 260 keV. The dashed line corresponds to the expected 
resonant part, the dotted line to the fitted non-resonant part, 
and the solid line is their sum including background. 



cascade transitions. The branching results are given in 
Table [I] they are in good agreement with previous work 
|l5j . except for the ground state transition. Our strength 
value of ui^y = 13.5 ± 0.4 (statistical) ±0.8(systematical) 
meV, which was determined absolutely, is also in good 
agreement with previous work 00] , ui^j = 14 ± 1 meV. 
The fitted total width of T = 0.99±0.03 keV is consistent 
with previous work ( 4]: V = 1.2 ± 0.2 keV). 

All transitions show isotropy at the En — 259 keV, 
J v = 1 /2 + resonance and no forward-backward asymme- 
try outside the resonance. Consequently, an angle- 
integrated 70-ray yield may be derived directly from the 
yield collected by the Ge detector at 9 1 = 55°. Finally, 
the summing due to the large solid angle of the detec- 
tor and the cascade coincident events in the transition to 
the ground state was corrected for by studying all other 
capture transitions. The absolute scale of the cross sec- 
tion deduced from the 7-ray line shape was obtained by 
normalizing to the thick target yield in the 259 keV res- 
onance 0] whose strength ^7 was determined absolutely 
(see above). 

The capture data from LUNA and previous work 4] for 
the transition to the 6.79 MeV state are shown in figure|3] 
there is an excellent agreement between both data sets in 
the region of overlap. As shown from previous work 00 , 
the S-factor has an external non-resonant contribution 
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FIG. 4: Full-energy peak for the 6791 — > keV transition 
observed with three Ge detectors positioned at 0° (open cir- 
cles), 90° (filled squares), and 125° (open triangles). The in- 
sert shows the measured energy shift relative to the expected 
full Doppler shift, where the solid line corresponds to the at- 
tenuation factor F(t) = 0.99. 



and a contribution from the Er = 259 keV resonance, 
where the latter contribution was derived from W7 and 
T values of the present work. The R-matrix fit 0, El 
included both contributions. The resulting ANC value of 
C = 5.0±0.1 fm _1//2 agrees with indirect determinations 
([13: C = 5.2 ± 0.7 fm" 1 ^; 0: C = 4.6 ± 0.5 fm" 1 / 2 ) 
and leads to the fitted curve shown in figure For large 
values of the R-matrix radius a the fit deviates from the 
data for energies above E = 1 MeV (see Fig. |5J). We 
have adopted the results for a = 5.5 fm leading to the 
extrapolated value S6.79 = 1.35 ±0.05 (statistical) ±0.08 
(systematical) keV b. This value is about 20 % lower 
than the R-matrix fit |f| of the data from |4j alone. 

The data from LUNA and previous work 0] for the 
ground state capture are shown in figure |SJ We have 
corrected the previous data for summing effects, which 
are at most 10 % above E = 500 keV, except near the 
destructive interference structure of the Er = 987 keV 
resonance. The data were fitted including the 3/2 + sub- 
threshold state, the 1/2+ , 259 keV, the 3/2+, 987 keV 
and the 3/2 + , 2187 keV resonances as well as a back- 
ground pole located at 6 MeV. For the subthreshold state, 
we used the reduced width 7 2 obtained from the fit of 
the data for the 6.79 MeV transition. The fit parame- 
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FIG. 5: Transition to the 6.79 MeV state in O. The S-factor data from present work are represented by solid points, those 
of by open triangles. The R-matrix fit (solid line) was obtained for a — 5.5 fm, the dashed line for a = 6 fm and the dotted 
line for a — 5 fm. The data point in the 259 keV resonance is off scale. 
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FIG. 6: Astrophysical S(E)-factor curve for the ground state transition in 14 N(p,7) 15 0. Filled-in data points are the results 
from LUNA, while the open data points are from previous work 4] correct for summing effects. The solid line corresponds to 
the R-matrix fit for a — 5.5 fm, the dashed line for a = 6 fm and the dotted line for a = 5 fm. 
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ters were the T 7 of the subthreshold state, the 987 keV 
and 2187 keV resonances, the T p and T 7 of the back- 
ground pole. For the external contribution we used the 
ANC of ^3 as a starting value. The results are shown 
in Fig. for a = 5 fin and C gs = 5.8 fin" 1 / 2 (dotted 
line), a = 5.5 fm and C gs = 7.3 haT 1 / 2 (solide line), 
and a — 6.5 fm and C gs = 8.8 fm -1 / 2 (dashed line). 
The corresponding T 7 values for the subthreshold state 
and extrapolated S-factor values are listed in Table ITU It 
can be noted that our deduced value T 7 = 0.8 ± 0.4 eV 
is in good agreement with the value from a life time 
measurement by [6] T 7 = 0.41+° 'fg eV as well as with 
T 7 = 0.95lo 95 e V, the value from coulomb excitation 
work 0. All T 7 values from Table |ITJ are consistent 
within uncertainties with our observation of a nearly full 
Doppler shift for the transition of the 6.79 MeV state 
(see Fig. 0}. The best x 2 for the fit to the ground state 
transition is obtained for a = 6 fm (x 2 = 1.27), however 
the differences are small (for a = 5 fm, x 2 = 1.35; con- 
tributions to x 2 near the narrow resonances have been 
omitted in the x 2 values given above). The best over- 
all agreement with all available data 0, 0, 0, El is 
obtained for a = 5.5 fm, where our fitted ground state 
ANC is C = 7.3 fm~ 1//2 , that can be compared with the 
result of [TjJ which gives C = 7.3 ±0.4 fm -1 / 2 , after con- 
version to the coupling scheme used in the present work. 
For the total S-factor a contribution from the transition 
to the 6.17 MeV state of S 6 .n{°) = 0.06 keV b from @ 
has been added to obtain an average (a = 5 to 6 fm) of 
S to t(0) = 1.7 ±0.1 (statistical) ±0.2(systematical) keV 
b, which can be compared with 1.77 ±0.20 keV b [5j and 
1.70 ± 0.22 keV b [17| from previous analyses. 
In summary, the present work improves the experimen- 



tal information concerning the ground state transition in 
14 N(p, 7) 15 0. The previous data 3 corrected by sum- 
ming effects discussed above are now in good agreement 
with the present work. The large ambiguity for the ex- 
trapolation of the ground state transition has been re- 
duced. In the energy region between E = 170 to 300 keV 
the S(E)-factor is dominated by the 259 keV resonance. 
Only at energies of about 100 keV and below the contri- 
butions of the subthreshold resonance and non-resonant 
mechanisms should become clearly observable. Using 
large volume Ge detectors it is impossible to obtain ad- 
ditional information on S SS (E) at lower energies (than 
current work) with acceptable uncertainty due to the 
sizable summing correction. An improved information 
on St t(E) can be achieved possibly using a 47r-summing 
crystal to observe the total S(E)-factor. Such an experi- 
ment is presently underway at the LUNA facility. 

After submission of the present work several relevant 
publications appeared: 

(i) The authors of reanalyzed the data of Q based 
on conclusions from their analyzing power data. This 
led to a different extrapolation for the transition to the 
6.79 MeV state (about 10 to 20 % lower than the analysis 
of [5| and 11 % higher than present work) and for the 
transition to the 6.17 MeV state (about a factor 3 higher 
than |5j). However, the primary aim of the present letter 
was to clarify the strength of the ground state transition. 

(ii) The analysis in [2(1 0, ^ were partly based on 
the results of the present work and led to the conclusion 
that the solar neutrino flux from CNO cycle is reduced by 
a factor 2 and the age of the globular clusters is increased 
by about 1 billion years. 
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